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(57) Abstract: The invention relates to a method for pre-qualifying subscriber connection lines for broadband services. According 
to said method, a mathematical characterisation of the physical behaviour of the subscriber connection line is obtained by applying a 
signal to the subscriber connection line to be tested and the output signal is then observed. Conclusions concerning relevant parame- 
ters of the subscriber connection line are inferred from said mathematical characterisation with the help of classical physical models. 
A combined time-frequency representation, the Weil transformation according to "A. Weil, Basic Number Theory, Grundlehren der 
mathematischen Wissenschaften, Bd. 144, Springer 1985", is used to obtain a statistically solid estimation of the aforementioned 
mathematical characterisation. 

(57) Zusammenfassung: Die Erfindung betrifft ein Verfahren zur Vorqualifizierung von Teilnehmeranschlussleitungen fur Breit- 
banddienste mittels dem nach Anlegen eines Signals an die zu testende Teilnehmeranschlussleitung durch Beobachtung des Aus- 
gangssignales eine mathematische Charakterisierung des physikalischen Vernal tens der Teilnehmeranschlussleitung erreicht wird. 
Ausgehend von dieser mathematischen Charakterisierung wird dann mit Hilfe von klassischen physikalischen Modellen auf rele- 
vante Parameter der Teilnehmeranschlussleitung zuriickgeschlossen. Zur statistisch robusten Schatzung der genannten mathema- 
tischen Charakterisierung wird eine kombinierte Zeit-Frequenz-Darstellung verwendet, die Weil-Transformation gemaB "A. Weil, 
Basic Number Theory, Grundlehren der mathematischen Wissenschaften, Bd. 144, Springer 1985". 
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VERFAHREN ZUM TESTEN VON TEILNEHMERANSCHLUSSLEITUNGEN FUR BRAEITBANDDIENSTE 

Verfahren zur Vorqualif izierung von Teilnehmeranschluss- 
leitungen fur Breitbanddienste 

Die Erfindung betrifft ein Verfahren zur Vorqualif izierung 
von Teilnehmeranschlussleitungen ftir Breitbanddienste. 

Die breitbandige Datenubertragung uber die Kupf erdoppelader 
der Teilnehmeranschlussleitungen des klassischen 
Telef onnetzes, insbesondere ADSL (asymmetrical digital 
subscriber line) , gewinnt zunehmend an Bedeutung. 
Gemeinsames Kennzeichen aller standardisierten xDSL-Dienste 
(ADSL, SDSL, SHDSL, VDSL) ist die Ausdehnung des genutzten 
Frequenzbereichs weg von der Sprachbandbreite (4Khz) bis 
hinein in Bereiche von uber 1 MHz . 

Da das Telefonnetz nicht ftir diesen Frequenzbereich . 
konzipiert ist, basieren diese Dienste auf Ratenadaptivitat 
der Ubertragungsverf ahren, d.h. die erzielbare Nutzdatenrate 
hangt jeweils von den technisch-physikalischen Gegebenheiten 
des individuellen Teilnehmeranschlusses ab. 

Haupteinf lussgrossen sind dabei die Leitungslange, der 
Leitungstyp und gegebenenf alls elektromagnetisch eingestreute 
Storungen; aber auch Stuckelungen, nicht- fachgerechte 
Hausverkabelungen, defekte Endgerate, offene Stichleitungen 
und Nebensprechen, wenn in einem Kabelbundel zahlreiche xDSL- 
Dienste auf geschaltet sind, verringern die moglichen 
Datenraten. 

Die nicht vorhersagbare Leistung (Datenrate) der 
Breitbanddienste bedeutet fur die Anbieter dieser Dienste ein 
hohes Erstinstallationsrisiko, so z.B. wenn der Anbieter bei 
der Installation eines derartigen Anschlusses finanzielle 
Aufwendungen auf sich nimmt, die er uber 
Internetnutzungsgebiihren zu erwirtschaf ten 
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beabsichtigt (Gratis-PC mit ADSL -Modem ) . Aber auch bei weniger 
riskanten Geschaf tsmodellen kann eine Abschatzung der 
moglichen Datenrate eines Breitbanddienstes jedenfalls 
Servicekosten reduzieren bzw. die generelle Kundenakzeptanz 
5 der XDSL-Technologie verbessern. 

Bekannte Testverf ahren fur Teilnehmeranschlussleitungen wie 
die sogenannte Zeit-Bereichsref lektometrie (kurz TDR fur 
time-domain ref lectometry) zum einseitigen, messtechnischen 
10 Test von Nachrichtenkabeln, wurden fur die Sprachbandbreite 
entwickelt und sind daher nur bedingt zur Vorgualif izierung 
fur xDSL-Dienste geeignet. 

Beim TDR-Verfahren wird ein kurzer (und damit relativ 
15 breitbandiger) Spannungspuls an das Adernpaar angelegt und 
die Verzogerung des Echos gemessen. Bei bekannter 
Ausbreitungsgeschwindigkeit des Signales kann daraus direkt 
die Kabellange ermittelt werden. Durch periodische 
Wiederholung der Messung und Oszilloskopie konnen von gut 
2 0 ausgebildeten Fachleuten auch Stuckelungen, offene 

Stichleitungen und ahnliche, die Datenrate beeinf lussende 
Storquellen erkannt werden. Diese Methode erfordert den 
Einsatz von gut geschultem Personal und ist damit aufwendig 
und kostenintensiv. 

25 

Der Erfindung liegt daher die Aufgabe zugrunde, ein Verf ahren 
anzugeben, das mit geringem Aufwand die Vorqualif izierung 
einer Teilnehmeranschlussleitungen fur Breitbanddienste und 
damit die genaue Vorhersage der erzielbaren Datenrate eines 
30 Breitbanddienstes erlaubt. Wesentlich ist es dabei, dass es 
das Verf ahren erlaubt, die Teilnehmeranschlussleitungen 
einseitig vom Amt aus zu messen. 

Gelost wird diese Aufgabe mit dem erf indungsgemaSen Verfahren 
35 nach Anspruch 1. 
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Die Erfindung wird anhand eines Beispieles naher erlautert: 

Kern des erf indungsgemafien Verfahrens ist ein System- 
identif ikationsalgorithmus mittels dem nach Anlegen eines 
5 Signals an die zu testende Teilnehmeranschlussleitung durch 
Beobachtung des Ausgangssignales eine mathematische 
Charakterisierung des physikalischen Verhaltens der 
Teilnehmeranschlussleitung erreicht wird. 

10 Ausgehend von dieser mathematischen Charakterisierung wird 
dann mit Hilfe von klassischen physikalischen Modellen auf 
relevante Parameter der Teilnehmeranschlussleitung 
zuriic kges chl ossen. 

15 Zur statistisch robusten Schatzung der genannten 

mathematischen Charakterisierung wird eine kombinierte Zeit- 
Frequenz-Darstellung verwendet, die Weil-Transformation gem&£ 
„A. Weil, Basic Number Theory, Grundlehren der mathematischen 
Wissenschaf ten, Bd. 144, Springer 1985". 



20 



Im folgenden wird das erf indungsgemafie Verfahren im Detail 
beschrieben: 

Es wird ein zeitdiskretes Multitrager-Zuf allssignal als 
Testsignal ausgesendet. Dieses ist wie folgt aufgebaut: 



25 




nk 



s(n) = %^c kJ g(n-lN T )exp j2n 



M 



F J 



mit folgenden Variablen: 



30 



n 



der Zeitindex (korrespondiert zum Sendetakt) 



s{n) 



die Abtastwerte des Testsignals (im Sendetakt) 



35 



der Symbolindex (korrespondiert zum Symboltakt ) 
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J ..der Tragerindex (korrespondiert zum Tragerf reguenz- 
abstand) 

C 

kJ komplexwertige Zufallskoef f izienten 

die Symboldauer (gemessen im Sendesignaltakt ) 

M 

f die Lange der der Modulation zugrundeliegenden 

diskreten Fouriertrans formation = (#Modulationstrager+l) *2 -1 

&( n ^ die (reellwertigen) Abtastwerte des Sendepuls 

(im Sendetakt) 

Die stochastische Innovation des Sendesignals rtihrt 
ausschliesslich von den komplexwertigen Zufallskoef f izienten 

her. Folgende Statistik wird dabei angestrebt : 

E{c u c mn }=a 2 {k)S km S ln , 

(2) 

d.h. die Koeff izienten sind unkorreliert und ihre Varianz 
hangt ausschliesslich vom Tragerindex ab ( ** m ist hierbei das 
20 Kroneckersche Delta-Symbol, und C mn ist die komplex 

konjugiert Zahl von C ) . 

Eine weitere, in Zusammenhang mit Basisband- 
Multitragermodulation durchaus iibliche Bedingung ist die 
hermitesche Symmetrie der Koeff izienten: 

(3) 

Dadurch ist sichergestellt dass das Sendesignal reellwertig 
ist. Gig. (2) kann naherungsweise durch einen ublichen 
3 0 Pseudozuf all s zahlengenerator (wie z.B. in W.H. Press, B.P. 

Flannery, S.A. Teukolsky, W.T. Vetterling, ^Numerical Recipes 
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in C w , Cambridge University Press, Cambridge 1993. 
und entsprechende Gewichtung gewonnen werden. 

Sofern die Periode des Pseudozuf allsgenerator grosser ist als 
5 die Messperiode, kann man das Sendesignal als Realisierung 
eines zyklostationaren Zuf allsprozesses betrachten. D.h. im 
Gegensatz zur klassischen Zeitbereichsref lektometrie ist das 
Signal nicht strikt periodisch (nur in seinen statistischen 
Eigenschaf ten gibt es Periodizitat ) . 

10 

Andererseits ist das Testsignal auch kein Spezialfall 
ublicher Multitragersignale (wie sie in der Trainigsphase von 
ADSL-Modems verwendet werden) die entweder einen zyklischen 
Prafix bzw. orthogonale Pulsformen gemSS S.D. Sandberg and 
15 M.A. Tzannes, ^Overlapped discrete multitone modulation for 
high speed copper wire", IEEE Journal Selected Areaes of 
Communications, Vol.13, No. 9, Dec. 1995. verwenden. 

Charakteristischer Bestandteil der Erfindung ist die 
2 0 Verwendung einer relativ langen (3-5 mal die Symbollange 
uberspannenden) gut f requenzlokalisierten 
Fensterfunktion, zusammen mit der Bedingung 

M F = KN T , 

(5) 

25 sodass letztendlich ein nichtorthogonales , ubervollstandiges 
Funktionensystem entsteht. 

Dies ist unterschiedlich zu dem Ansatz von W. Kozek und A. F . 
Molisch, „Nonorthogonal pulseshapes for multicarrier 
30 communications in doubly dispersive channels", Oct. 1998, 
IEEE J. Select. Areas Comm., Vol.16, No. 8, pp . 1579-1589 . 
Bei dem ein nichtorthogonales, aber unvollstandiges 
Funktionensystem fur die Multitragerubertragung vorgeschlagen 
wird. ) 

35 
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Insbesondere bedeutet es, dass die Redundanz des Sendesignals 
gewis sennas sen iiber Zeit- und Frequenzrichtung verschmiert 
ist. Solche Funktionensysteme sind lediglich bekannt aus der 
Signalanalyse, die Anwendung zur Syn these eines 
teilredundanten Pseudozuf allssignals allerdings neu und 
erf inderisch. Dieses erf indungsgemafie Sendesignal bietet 
folgende Vorteile: 

Die Messbandbreite kann flexibel gestaltet werden, d.h. es 
konnen zum Beispiel die POTS-Bander (also der fiir klassische 
Telefonie benotigte Frequenzbereich ausgeblendet werden 
(notwendig fur die Integration im 

Breitbandanschlussmultiplexer bei ADSL) . Dies ist ein klarer 
Vorteil gegenuber W.W. Jones, „ Sequence Time Domain 
Reflectometry (STDR) For Digital Subscriber Line Provisioning 
and Diagnostics" , White Paper, Mindspeed Technologies TM, wo 
ein sich von der Gleichspannung beginnender Frequenzbereich 
verwendet wird. 



Die verfugbare Messbandbreite kann wesentlich besser 
ausgenutzt werden. Die Empf indlichkeit gegenuber 
Pulsstorungen ist geringer. Bei nicht streng synchroner 
Sende- und Empf angssignalabtastung kann der Frequenzof f set 
explizit kompensiert werden. 

Durch Phasenmittelung kann man den zyklostationaren Prozess 
auf einen stationaren Prozess abbilden. Das 

S ( f) 

Leistungsdichtespektrum sKJ } des phasengemittelten 
Testsignals hangt von der Lei stungsvertei lung auf den Tragern 

ab (durch G ^ festgelegt) und der Wahl des Pulses . 
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V 



f -f 



J 



(6) 

wobei T die Abtastperiode des Sendesignals ist und ^^^das 
Spektrum des Sendepuls. Durch entsprechende Auswahl von 

flr (^) konnen die zulassigen Masken fur das Sendespektrum zum 
Beispiel Ausblendung der POTS-Bander erreicht werden. Die 

Wahl des Pulses S( n ) ± st wen iger kritisch, im einfachsten 
Fall wahlt man eine klassische, gut f requenzlokalisierte 
Fens terfunkt ion wie zum Beispiel Hamming oder Kaiser, siehe 
W. Kozek und A.F. Molisch, „Nonorthogonal pulseshapes for 
multicarrier communications in doubly dispersive channels", 
Oct. 1998, IEEE J. Select. Areas Comm., Vol.16, No. 8, 
pp. 1579-1589. 
Seite 126. 

Im weiteren werden Auschnitte der Lange N aus dem 
Sendesignal wie folgt betrachtet: 

* m in) = s(n - mKN T )Z [W (n - mKN T ) 

(7) 



hierbei steht 
und Lange N 



fur einen Rechteckpuls mit Amplitude 1 



Dae Weil -Trans formation eines Vektors 

x(n) 



der Lange 



N = KN. 



ist gegeben durch 



25 



W x (n, k) = -j= y Z X (n +IN T ) exp 




(8) 
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sie steht in umkehrbar eindeutiger Beziehung zum Signal ( A. 
Weil, Basic Number Theory, Grundlehren der mathematischen 
Wissenschaften, Bd. 144, Springer 1985) . Physikalisch kann 

hier m als Zeitindex und ^ als Frequenzindex 
interpretiert werden. Grundprinzip des Identif ikations- 
verfahrens ist nun iiber eine Kreuzkorrelationsanalyse im 
Weil-Bereich zu einer Schatzung der Weil-transf ormierten 
Echo-Impulsantwort zu kommen. Die Kreuzkorrelationsfunktion 
ist hier definiert als: 



C y§s (m, k) = E^V y (n + m, k)W s (rc, k)\ . 



'y 

(9) 

Vorteilhaf terweise kann die Weil-Transf ormierte eines 
Sendesignalblockes durch zweidimensionale, diskrete 
Fouriertrans formation der Sendekoef f izienten und 
Multiplikation mit der Weiltransf ormierten des 
Sendesignalf ensters gebildet werden 

w 



N T i+K-1 

2 2 

m = 0 I = i 



c kJ exp 



- j2x 



mk 



v 



N 



+ 



nl 

IF 



J 



) 



(10) 



Weiters werden Empf angssignalkoef f izienten durch ein 
Skalarprodukt mit dem modulierten Empf angssignalf enster wie 
folgt gebildet: 

(l+l)N T _ s 

<*k,l= E ytoTO* ~ W T ) exp - j2n 
n=lNj \ 



nk 
Mr 



Daraus kann dann analog die Weiltransf ormierte des 
25 Empfangssignals y(n) analog zu (10) gebildet werden: 



WO 03/088574 



» 



PCT/EP03/03930 



W 



yZ[i,i+KN r = W y (m,n)- 



N T i+K-l 

E E d kje*P 
m = 0 / = z 



- j2n 



mk nl 

■ + 



v 



Unter der Annahme der Ergodizitat des Eingangs signals (nur 
gerechtfertigt sofern die Periode der beteiligten 
Pseudozufallssignale grSsser als die Messzeit bleibt) , kann 
man die Erwartungswertbildung durch Zeit/Frequenzmittel- 
wertbildungen ersetzen: 

i N M 

Cy, s (m,k)=— — — £ W y (n t k)W Sm in + m,k), 

(11) 

wobei ^^den in Gig. (7) beschrieben Sendesignalblock 
beschreibt . 



Das vorgeschlagene Verfahren zur Systemidentif ikation eignet 
sich natiirlich auch zur adaptiven Echounterdriickung (echo 
cancelling) , mit den oben beschriebenen Vorteilen der 
Moglichkeit von Ausblendung bestimmter Frequenzbander , wie 
beispielsweise in einer der folgenden Literaturstellen 
beschrieben : 

K. Townsend et al, „Apparatus and method for echo characteri- 
zation of a communication channel US Patent Nr. 5,577,116, 
Nov. 1996. 



J. P. Agrawal et al, „Adaptive echo cancellation and equali- 
zation system signal processor and method therefor", US Pat- 
ent Nr. 4,760,596, July 1988. 

S. Wu et al, „Efficient echo cancellation for DMT MDSL", US 
Patent Nr. 6,072,782, June 2 000. 
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M. Ho et al, „Method and apparatus for echo cancellation 
with discrete multitone modulation" , US Patent Nr. 5,317,596, 
May 1994. 

Die beschriebenen Verfahren zur Systemidentif ikation (als 
Bestandteil eines adaptiven Echo-Cancellers) unterscheiden 
sich in mehreren Gesichtspunkten vom erf indungsgemaSen 
Verfahren, zumeist werden Abwandlungen des sogenannten LMS- 
Algorithmus verwendet oder das klassische Kalmanf ilter . 

Nach der vorliegenden Anwendung zur Vorqualif izierung von 
Teilnehmeranschlussleitungen fur Breitbanddienste geschieht 
die Auswertung der Weiltransformierten der interessierenden 
Impulsantwort wie folgt: 

Zwar kann bereits ein geubter Betrachter aus Rohdaten (z.B. 
der Echoimpulsantwort) einige Riickschlusse auf die Lange 
einer Leitung bzw. die Art des Absehlusses ziehen, fur den 
Einsatz der vorliegenden Erfindung ist aber eine 
automatisierte Auswertung vorgesehen, deren Endergebnis die 
Vorhersage von xDSL- Performance ist ( >x reach prediction^). 

Etwaige additive Storungen die fur die Ermittlung der xDSL- 
Performance unter Umstanden relevant sind, konnen direkt aus 
den ohnehin schon ermittelten Emp fangs signalkoef f izienten 
durch Erwartungswertbildung uber den Zeitindex gebildet 
werden : 



S = Er < 

noise , k 



. 2 
d l,k 



Die dabei gewahlte Vorgangsweise besteht nun darin, ausgehend 
von gemessenen Rohdaten (komplexe Zahlen) eine multi- 
dimensionale Tabelle anzulegen, die durch relevante 
physikalische Grossen parametrisiert ist (z.B. Lei tungs lange 
und Abschluss) . 
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Im folgenden wird die Vorgangsweise fur die Parameter 
Leitungslange und Abschluss beschrieben, die Auswertung 
konnte aber konzeptuell analog auch auf andere 
Zwischenparameter abzielen. 



W (v a) 

Die Rohdaten h\r*is werden mit als Tabellenvektoren 

T (k m) (v a) 

dargestellten Werten von gemessenen 
Referenzleitungen verglichen und es werden jene Indizes 
ausgewahlt, die die kleinste absolute Abweichung (im Sinne 
komplexer Zahlen aufweisen) : 



(C,L)= argminE ^T^\p,q)-W h {p, q ) 

y y kj q P 



Aus diesem Vergleich kann damit direkt auf die physikalischen 
Eigenschaf ten der zu messenden TeilnehmeranschluSleitung 
geschlossen werden. 

Hierbei entspricht der Index k der beispielhaf ten 
Implementation folgend einem Widerstandswert des Abschlusses 
und der Index 1 einer Leitungslange. Urn die Auflosung des 
Messverfahrens zu erhohen, erweist es sich als vorteilhaft 
nach dieser Optimierung eine komplexe Interpolation 
anzuschliessen. Diese kann mit den liblichen Verfahren zur 
multidimensionalen Interpolation durchgefuhrt werden, wie sie 
z.B. in W.H. Press, B.P. Flannery, S.A. Teukolsky, W.T. 
Vetterling, ^Numerical Recipes in C M , Cambridge University 
Press, Cambridge 1993 beschrieben sind. 



Die letztendlich interessierenden xDSL Parameter (zu 
erwartende Datenubertragungsrate ) konnen dann durch eine 
einfache Tabellensuche gewonnen werden. 
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Eine Kalibrierung des Messsystems ist dann zweckmaSig, wenn 
z.B. bei langen Leitungslangen das interne Echo der 
Messanordnung typischerweise wesentlich starker sein wird als 
5 das charakteristische Echo vom Ende das Kabels. Beim 

Ausfuhrungsbeispiel ist eine automatische Kalibrierung durch 
Anschaltung digital steuerbarer Leitungssimulatoren bzw. 
Widerstandsdekaden durchf uhrbar . 



10 
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Patentanspruche 



10 



15 
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1. Die Erfindung betrifft ein Verfahren zur Vorqualif izierung 
von Teilnehmeranschlussleitungen fur Breitbanddienste mit 
folgenden Verf ahrensschritten: 

- an einen Testpunkt der TeilnehmeranschluSleitung wird ein 

M N f flk 

gemafi ^) = EE c u ^ ,2 "^) eX P J 2?t ' 



*=o /=o 



gebildetes 



zeitdiskretes Multitrager- Sendesignal angelegt, 

die Echoimpulsantwort y(n, der TeilnehmeranschluSleitung 

wird an dem Testpunkt gemessen, 

aus der Echoimpulsantwort werden zweite komplexwertige 
Zuf allskoef f izienten d kfl gemaS 



« = (/ + i)W r 



d kj = X y(n)y(n - IN T ) exp 
« = (/-i)^ r 



- j2n 



nk 



M 



J 



ermittelt, 

der empirische Schatzwert der Kreuzkorrelationsfunktion 

WAp,q) . 

der aus den Indizes c kil md d kll durch 
zweidimensionale diskrete Fouriertrans formation gebildeten 
Signale 



N T i+K-1 

Cm]n =2 X C k,l ex P 
m = 0 l = i 



r 

- j27Z 



V 



mk nl 
+ — 
K 



N 



N j- i + K - I 

D m]n = X X d k,l 
m — 0 / = i 



r 



exp 



r 



- j2n 



mk nl 
N T ~K 



J 



J 



J 
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gemaS (wobei 0</l<l ein Vergessensf aktor ist, der an die 
gesamte Mittelungslange der Messung entsprechend der 
Rechengenauigkeit des verwendeten Prozessors zu wahlen ist) 

K <|+» = (1 - X)K^ k + XC<t\ n , k W t (m + n,k)D^ n W- r (». 



wird ermittelt, 
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und analog zur Schatzung der Kreuzkorrelat ions funkt ion eine 
Schatzung des Leistungsdichtespektrums etwaiger Storsignale 
entsprechend 

5 durchgefuhrt , 

- der empirische Schatzwert der Kreuzkorrelationsf unktion 

W (p q) 

h\r>*is wird mit den gespeicherten Werten von gemessenen 

Referenzleitungen v/ 7 **// verglichen, und aus dem 

Vergleich werden die physikalischen Parameter der 
10 TeilnehmeranschluSleitung ermittelt. 
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